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Entangled qubit pairs can serve as a quantum memory or as a resource for quantum communi- 
cation. The utility of such pairs is measured by how long they take to disentangle or decohere. 
To answer the question of whether qubit-qubit interactions can prolong entanglement, we calculate 
the dissipative dynamics of a pair of qubits coupled via the exchange interaction in the presence 
of random telegraph noise and 1/f noise. We show that for maximally entangled (Bell) states, the 
exchange interaction generally suppresses decoherence and disentanglement. This suppression is 
more apparent for random telegraph noise if the noise is non-Markovian, whereas for 1/f noise the 
exchange interaction should be comparable in magnitude to strongest noise source. The entangled 
singlet-triplet superposition state of 2 qubits {'ip± Bell state) can be protected by the interaction, 
while for the triplet-triplet state {4>± Bell state) , it is less effective. Thus the former is more suitable 
for encoding quantum information. 

PACS numbers: 03.67.-Pp, 03.65.-Yz, 05.40.-a 



I. INTRODUCTION 

Much theoretical and experimental effort has been di- 
rected towards studying the viability of quantum infor- 
mation processing (QIP) in recent years due to a se- 
ries of remarkable QIP algorithms [1-3^. Initial concerns 
about quantum coherence being too fragile to be use- 
ful have been partially dispelled with the discovery of 
quantum error-correcting codes [IHZj, quantum thresh- 
old theorems [THTU]. decoherence- free subspaces [TIHT5] 
and dynamical decoupling by using optimized pulsed se- 
quences [HHlT]. 

Qubit decoherence can be attributed to various sources 
and has been investigated using models such as spin- 
bath models [T8Vl20j . hyperfine interaction models [2T1 - 
[25] and phonon induced decoherence pilES] . Another 
common source of decoherence in solid state devices are 
two level systems (TLS) that generate random telegraph 
noise (RTN) , which with a wide distribution of switching 
rates can give rise to 1/f noise [171 HB]. In a number 
of recent semiconductor quantum dot(QD) experiments, 
RTN is observed when the potential in the dot lines up 
with the electrochemical potential in the reservoir, caus- 
ing electrons to randomly tunnel back and forth between 
the dot and the reservoir p5^31j . This puts limitations 
on the performance of a QD qubit, as such a random 
telegraphic current can modulate the QD's orbital wave- 
function which can create magnetic noise via spin-orbit 
coupling mechanisms. This type of noise is also known to 
be important in superconducting qubits [321 , and could 
affect the performance of other types of qubits as well 

[3HMj- 

This suggests that ways be invented to prolong quan- 
tum coherence and entanglement in qubit pairs in the 
presence of such noise sources. In this paper we present 
an extremely encouraging set of results - that the ex- 
change interaction between the qubits can be used to 
suppress decoherence as well as disentanglement due to 
RTN. This is a natural choice as the Heisenberg exchange 



interaction between the qubits is often used in any case 
to implement various gate operations such as controUed- 
NOT and SWAP gates [37H40|. Thus new circuit ele- 
ments are not required. Typical proposals to suppress 
decoherence from TLS that rely on spin echo techniques 
[STl [2H1 im 112] require additional resources and system 
monitoring. Our proposal can be used either as an alter- 
native or even in addition to pulsing. 

Our aim in this paper is to show that for the maxi- 
mally entangled Bell states, the interaction between the 
qubits can be used to suppress decoherence and disen- 
tanglement. Without any loss of generality, we first ana- 
lytically and numerically show this effect using a model 
where a single RTN source is coupled to only one of the 
two qubits. We then consider a model with two uncorre- 
lated RTN sources, with each of them coupled to a qubit 
and show that the interaction suppresses quantum dis- 
sipation in this case as well. Our analytical results sug- 
gests that a straightforward generalization can be made 
to the case of multiple uncorrelated RTN sources. We 
subsequently show that the qubit interaction suppresses 
decoherence even if the qubits are coupled to a large num- 
ber of uncorrelated fluctuators with a 1// noise power 
spectrum. These results are extremely important for 
QIP as much of its vaunted capabilities are due to the 
fact that unlike the classical bit, multiple qubits exhibit 
quantum entanglement which allows multiple states to 
be addressed simultaneously[43 . Alternatively, an effec- 
tive single qubit can be created from the exchange cou- 
pled singlet-triplet states, which will be less susceptible 
to RTN. 

The exchange interaction in our work is taken to be 
an externally controllable parameter. Testolin et al. |33] 
have proposed a model to show how the two level fluctua- 
tors itself affect the exchange interaction as a function of 
time, while Das ^45] has shown that the interactions can 
lead to periodic disentanglement and entanglement be- 
tween the qubits in contact with different environments. 

In general, non-Hermitian Hamiltonians are often used 
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to describe decay processes in open quantum systems 
[46H50] . Our calculations are done using a recently de- 
veloped quasi-Hamiltonian formalism which is suitable 
for describing the non-unitary temporal evolution of a 
quantum system acted on by a classical stochastic pro- 
cess [Ml EU [52] ■ Similar approaches have been used in 
Refs. |531 154) . In many instances, Quantum and classical 
noise models can yield the same solution for decoherence 
if the noise power spectrums have the same line-shape. 
For example, in [55j [56] the analytical expression ob- 
tained for a qubit's decoherence in the presence of a ther- 
mal bosonic reservior in a lossy cavity with Lorentzian 
broadening is exactly the same as what we obtain for clas- 
sical pure dephasing noise. More recently, Saira et a^[53] 
have shown that within the Born approximation, many 
fully quantum mechanical noise models can be exactly 
mapped onto classical stochastic noise models, including 
those for RTN. 

This paper is organized as follows. First we describe 
our single fluctuator model and the quasi-Hamiltonian 
method in secjll] Then in sec |III[ we give results and the 
exact Bloch vector solutions for some important cases. 
Since for QIP applications the entanglement dynamics at 
immediate times is often the most important, in sec |IV| 
we obtain analytical results for the short time behav- 
ior for the Bell states for any arbitrary set of noise pa- 
rameters. Our results are then extended to the case of 
two uncorrelated fluctuators in secO This treatment of 
two fluctuators is then subsequently extended to treat an 
arbitrary number of uncorrelated fluctuators in sec |VI[ 
where we calculate the temporal dynamics of two inter- 
acting qubits in the presence of eight RTSs with a 1// 
distribution. Finally we present our conclusions and im- 



plications for qubit design in sec VII 



II. MODEL AND METHOD 

It is well known that in various materials and systems 
sudden step like transitions occur at random intervals 
of time between two or more discrete voltage levels [57]. 
In the case of semiconductors, these random telegraphic 
signals (RTS) are often attributed to trapping and release 
of charge carriers by defect sites. If the fuctuators are 
statistically independent, then the RTS can be expressed 
as a sum of the contributions from individual fluctuations 



(1) 



where, Si{t) is a RTN sequence that switches between ±1 
at random intervals of time, gi is the noise vector for the 
i*^ fluctuator and |g| = g is the noise strength. Their au- 
tocorrelation function is {si{ti)sj(t2)) (x exp(— 27^1^1 — 
t2\)Sij and 7i is the switching rate of the fluctua- 
tor. Each individual fluctuator has a Lorentzian power 
spectrum and a broad distribution of 7 results in a 1// 
noise power spectrum[57 . Analytical solutions are how- 
ever more tractable for few RTN sources, hence in this 



paper we first primarily focus on the dissipative effects 
of a single fluctuator. We then discuss the two fluctua- 
tor case, the results of which are then extended to treat 
a larger number of uncorrelated fluctuators with a 1// 
noise power spectrum. 

We use the following Hamiltonian to describe the two 
qubit quantum system. 



H 



Hn + Hn 



Hi, 



B • [Si + S2] + s{t)g ■ Si + JSi • S2 



(2) 
(3) 



where. Si — I(8)cri and S2 — <T2(X)I respectively repre- 
sents qubits one and two, cr is the triad of Pauli matri- 
ces, B is the steady magnetic field chosen to be in the 
z direction for all our calculations, and J is the Heisen- 
berg interaction strength. Note that in this model the 
RTN is only coupled to one qubit. Earlier calculations 
on non-interacting qubits suggest that this is sufficient 
to describe all the qualitative effects (55]. The angle 9, 
between the noise vector g and magnetic field B is called 
the working point of the qubit. The Hamiltonian H is 
written for a given realization s (t) of the noise. Physical 
quantities are calculated by averaging over all sequences. 

Since we are primarily interested in disentanglement 
of the qubits, it is convenient to rewrite the two qubit 
Hamiltonian using the maximally entangled Bell states 
V-i = (|10) ± |01)) /^/2 and (/i± = (|00) ± |11)) /^/2 as 
our basis states. Eq|3]in the 4., ?/;_, -0+] basis is 



^0 ^xy 



(4) 



gxs{t)I 



where = JI - [B^ + gzs{t)](Tx, H^y = 
igys{t)<Tx and = -[JI -|- 2Jaz + 5zs(i)cra;]- 

In general, density matrices are suitable for treating 
open quantum systems where one has a statistical mix- 
ture of pure states. For two qubits, the time dependent 
density matrix can be expressed as 



Pit) = 



1 



(»J)5^(0,0) 



nij[t)(7i (g) a J 



(5) 



the coherent time evolution of which is governed by the 
von Neumann equation. Here / is the 4x4 identity 
matrix, n(i) is the fifteen-component generalized Bloch 
vector. We have used a notation, where two indices i and 
j (where i,j = Q,x,y,z and 7^ (0,0)) are used to 

denote each component of the generalized Bloch vector. 
Here, ai^aj = are the generators of SU (4) and ctq = I 
is the 2x2 identity matrix. |n| can be thought of as a 
measure of purity. For instance, |n| = is the completely 
mixed state. 

In the recently-developed quasi-Hamiltonian formal- 
ism, the dissipative temporal dynamics of an open quan- 
tum system under the influence of classical noise is calcu- 
lated by transforming a random time-dependent Hamil- 
tonian into a time-independent non-Hermitian Hamilto- 
nian [3ni[5IllSl]- This is done as follows. For a given noise 
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realization the density matrix is unitarily time evolved as 
p{At) = U ■ p(0)•C/^ where U = exp{-iH{t)At/h). Sub- 
stituting these in Eqisland using the identity Tr {XiXj) — 
A5ij [58| , one obtains the following temporal transfer ma- 
trix equation 



applied to describe dissipation, while in the latter case 
the methods such as the quasi-Hamiltonian method is 
required. 



III. RESULTS AND DISCUSSION 



N 



n(i)- [] T„, -11(0). 



(6) 



where, N = t/ At and the transfer matrix elements are 
given by = Tr [UXiWX-j] /2. As the RTN is modeled 
as a classical stochastic process, its dynamics is governed 
by the master equation, W(i) = VW(f) [59], where V 
is a matrix of transition rates (such that the sum of each 
of its columns is zero) and W is the flipping probability 
matrix for the TLS. If the average occupation of the two 
states is the same (for unbiased fluctuators) , then 



W(<) = 



1 



(7) 



Here, 7 is the switching rate for the TLF. The combined 
temporal dynamics of the quantum and classical TLS, 
averaged over all noise sequences, can be described by 



nit) = (/|r^ |z)n(0) = (/I exp(-zF,t)|z)n(0) (8) 

where, \i) and |/) are the initial and final state vectors 
for the TLS that satisfy W|i(/)) = \i{f)). For an un- 
biased TLS (i.e. with equal occupation probabilities), 
\i) = I/) = [l,l]/y2. Here = rArrAr_i...ri and in 
the small time approximation, any of the matrices 



W0T 



(1 - 7At)T,=+i 7AtT,=+i 
7AtT,=_i (1 - 7At)T,=_i 



(9) 



Where denotes a Hadamard product and T is a square 
matrix, each of whose columns consists of the transfer 
matrices [T^^+i, Ts=-i]. T is now a 30 x 30 matri- 
ces in the combined generalized Bloch vector and TLS 
spaces. Hq = lim z(r — I)/At is a time-independent 

non-Hcrmitian quasi-Hamiltonian|52| . As Hg is not her- 
mitian, the time-evolution operator exp(—iHqt) is not 
unitary, which makes it suitable for treating dissipative 
processes in open quantum systems. 

There are thus two relevant time scales that determine 
the noise characteristics of the system. The correlation 
time of the environment (which is proportional to 7"^) 
and the time period for the noise induced Rabi-like os- 
cillations of the qubit ( which is proportional to g^^). 
If the noise correlation time of the environment is much 
less than the time period of the noise induced Rabi-like 
oscillation {g << 7, ), we say that the noise is Marko- 
vian. Whereas if 5 >> 7 we say that the noise is non- 
Markovian. In the former case Redfield theory can be 



At the pure dephasing point (i.e., 6 — 0), H^y — 
in Eq|4] and the two qubit dynamics can be effectively 
decoupled into two single qubit problems, whose dynam- 
ics is governed by and H^. The effective one particle 
quasi-Hamiltonian for ^± can be extracted using the pro- 
cedure outlined in the previous section. 



q 



H 



qj 

CTo) ® Lf) 



(10) 

- gzCTz ®Lx]-i [2J(To L^] 



where, L^^y.z G 5*0(3) and Lq is the three dimensional 
identity matrix. In the absence of J, the exact solution 
for the non-zero component of the Bloch vector using 
Eq|8]is 



nz{t) 



cos 



-It 



(11) 



where Q = y 5^ — 72. Note that as one crosses over 
from the non-Markovian to the Markovian noise regime 
(7 > gz), the trigonometric functions in nz{t) become 
hyperbolic functions and the Bloch vector's oscillations 
between -0+ and ip- would then not be seen. When 
5^ << 7 this reduces to the usual Redfield form nz{t) w 
exp(-i/T2) where, T2 = 7/52. 

For the effect of J on the Bloch vector, we approximate 
the matrix exponential in EqjS] using the Zassenhaus ex- 
pansion jSni as follows 



iH 



f q.l p^^ q(J ' q.J i / 



90 e 



(12) 



This gives the following expression for the Bloch vector 
component, valid for either small J, g^ or at short times 





FIG. 1: (color online) Effective Bloch spheres and the 
ensemble averaged Bloch vector trajectory at the pure 
dephasing point shown for (a) xp± and for (b) (j)± with 
= 1. The free precession of Joz about 7/1+ is also shown 
in (a). 
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n'^{t) w C cos( Jg^t^) + ^ sin(rJt) sin(2Jt) sin(J5^t2) 



-7i 



(13) 

where C — cos{ftt) + ^ sin(r2i). Now, in case of the two 
qubit Bloch vector for the non zero components are 
n-xx = TT-yy = Uz and Uzz = — 1- If in addition the ini- 
tial state hes entirely in the ^p± subspace, then the entire 
history of the qubit pair can be visuahzed in the cor- 
responding effective Bloch sphere, with ip± at the poles 
(see Figjl]). 

Let us first consider the dynamics on the 'ip± subspace, 
taking ip^ as the initial state. In the absence of J in 
the strong coupling limit, the effective Bloch vector os- 
cillates between -0+ and "0- with the quantum coherence 
and entanglement dissipating in time (see Figs{l]-a and 
[2]-a). This oscillation between ip^ and V'- is due to the 
gzs(t)ax component of H^. In a given noise realization, 
s(t) causes rotations about the x-axis which switch ran- 
domly between the two orientations. However the en- 
semble averaged Bloch vector always travels in a straight 
line from pole to pole for ^± (as shown in Fig{T]-a) and 
eventually diminishes to the center, since the averaging 
restores the chiral symmetry. In Fig. [2] we show |n| , the 
magnitude of the Bloch vector and the concurrence as a 
function of time and J for the 0+ Bell state. For bipar- 
tite systems, the concurrence provides a measure of the 
entanglement between two qubits and is defined as [61j 



[V^i - - n/ks - ^^■4, 0] 



(14) 



where, ki, K2, K3 and K4 are the eigenvalues of p{(Jy <Si 
'^y)p*{'^v ®'^y) in decreasing order. |n| , as stated above, 
provides a measure of the purity. The two quantities tend 
to track each other but are not in one-to-one correspon- 
dence. The dissipative dynamics of the entangled qubits 
is shown in three different noise coupling regimes -in the 
strong coupling limit (non-Markovian noise, g > 7), in 
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FIG. 2: (color online) Bloch vector magnitude and 
concurrence as a function of time and Js for shown for 
(a) non-Markovian, (b) intermediate and (c) Markovian 
noise coupling regimes at S = 0, and _Bz = 1. Note that the 
Bloch vector amplitudes are offset by -1-1 for clarity. 7, 
g and J are in units of B. 
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FIG. 3: (color online) Same as Figjijbut at 6 = 7r/4 



the intermediate regime {g ~ 7) and in the weak coupling 
limit (Markovian noise g < -f). Note that the purity of 
the state as measured by \n\ is zero only at the center of 
the Bloch sphere whereas the concurrence is zero every- 
where on the straight line that passes through the center 
and connects |01) to |10) (i.e. the x-axis). 

As the exchange interaction is turned on, the Jaz com- 
ponent of Hj, causes the Bloch vector to precess about 
-0+ (see FigjUa). This effect competes with the effect of 
the noise, which is to drive the system toward the ori- 
gin of the effective Bloch sphere. Therefore the effective 
ensemble averaged Bloch vector tends to remain closer 
to 0+ with increasing J; this delays decoherence and 
the disentanglement. This behavior holds true even in 
the intermediate noise and weak noise coupling regimes 
as shown in Figsj2]-b and c. However in the Markovian 
limit this suppression of decoherence occurs only at larger 
J values, namely when J ~ 7. This is shown in Fig[2]-c. 
The overall dissipation is also much slower for Markovian 
noise. 

We next consider an equal mixture of dephasing and 
relaxational noise, a working point of 6* — 7r/4.. The re- 
sults are shown in Figj3j Here no analytic solution is pos- 
sible even for J = and we must diagonalize the quasi- 
Hamiltonian Hq numerically. Qualitatively, the temporal 
behavior of dissipative process is the same as that of the 
pure dephasing case for all three noise coupling regimes, 
but with overall longer decoherence times. 
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FIG. 4: (color online) Time period for the (a) envelope 
function decay and (b) oscillations as a function of J and t 
J is in units of B. 
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Clearly, the presence of the interaction J suppresses 
disentanglement and decoherence for the initial state. 
In the non-Markovian case, however, there are several 
time scales, one for the envelope function decay, which 
fits well to a form ~ exp (—t/Tenv) and one for the oscil- 
lation period Tosc- In Fig|4j we show T^nv and Tqsc as a 
function of 9 and J. 

In Figj4]-(a) it is seen that Te„^ increases monotonically 
as a function of 6. On the other hand, T^nv initially de- 
creases as a function of J, reaches a minimum and then 
increases monotonically. This non-monotonic behavior is 
easily understood physically. Slow dephasing noise (small 
9) tends to mix ip- in with the initial state moving 
the Bloch vector's trajectory "south" on through the ef- 
fective Bloch sphere. Since the noise is non-Markovian, 
this can actually cause oscillations between -0+ and V'-- 
On the other hand, J tends to move the trajectory "east" 
or "west" for a fixed realization of the noise, causing the 
oscillation to just miss the south pole. This results in 
the speeding up of the oscillations on averaging. For suf- 
ficiently large J, however, the trajectory tends to stay 
entirely in the northern hemisphere, which slows down 
the oscillations. Here we have the basic mechanism for 
exchange interaction induced suppression of decoherence. 
If the initial state is an eigenstate of the interaction and 
the noise connects this state to a state belonging to a 
different eigenstate with a different eigenvalue, then the 
interaction term causes the trajectory to undergo tight 
oscillations near the initial state. 

To show this, we look at how the the action of J de- 
pends on the initial state. Take 0+ as the initial state. 
Then the question is how this mixes with (/>_. The cru- 
cial difference is that and (p- both belong to the 
triplet manifold, i.e., the exchange interaction has the 
same eigenvalue for the two states. 

When 9 = 0, (pure dephasing), the qubit's dissipative 
dynamics is independent of J, as seen in (Eqj4]) . The 
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FIG. 5: (color online) Bloch vector magnitude and 
concurence for 0+ as a function of time and Js shown for 
(a)non-Markovian, (b)intermediate and (c)Markovian noise 
coupling regimes at 9 — 7v/6 and at = 10~^. The Bloch 
vector amphtudes are offset by -1-1 for clarity. 7, g and J are 
in units of B. 




lgl=0.1.7=0.1 

^ J=0.1 
\Vj=0 




(b) 



■ |n(t)| 
■C(t) 



(c) 



O 50 100 150 50 100 150 50 100 150 
Time (1/B) Time (1/B) Time (1/B) 

FIG. 6: (color online) Same as FigjHjbut &t 6 — tv/3. 



effective one particle quasi-Hamiltonian is therefore 

= I [7(a^ - ao) ® + {2B,ao + ga,) (g) L^] . (15) 
The exact solution of which gives 



ny{t) 



cos(m) + ^ sin(m) cos(2B^t)e-^* (16) 

n^{t) = - cos(m) -I- ^ sm{nt) sin(2S^i)e"''* (17) 

The ensemble averaged temporal trajectory of the Bloch 
vector at a large magnetic field is shown in Fig{T]-b. In 
the absence of a magnetic field, the Bloch vector only 
travels in a straight line from pole to pole. However, as 
in the case of '0±, the overall magnitude of the Bloch 
vector (hence the decoherence) for (f)± is independent of 
the magnetic field at 6* = as seen from EqjlT] The cor- 
responding non-zero two qubit Bloch vector components 



are, ri; 



xy 



and n^. 



^yy 



If 7^ 0, then Bz suppresses the decoherence. How- 
ever the effects of J are not visible if i?2 >> J, hence we 
set Bz = 10~^ for the next set of calculations. In Figsjs] 
and [6) the decoherence and disentanglement dynamics 
is shown for 0+ as a function of time and J, at various 
noise coupling regimes for 9 = tt/6 and 9 = tt/3 respec- 
tively. In the absence of J, the decoherence has the same 
qualitative behavior as that of Eq. [17] With the onset 
of J the decoherence is suppressed initially, however at 
longer times it exhibits a crossover behavior in the non- 
Markovian and intermediate noise coupling regimes. This 
crossover occurs later in time with increasing 9. Thus the 
picture is more complicated than in the case of ^± and 
the beneficial effect of the exchange interaction is weaker. 

The results presented in this paper using the quasi- 
Hamiltonian method have also been confirmed through 
numerical simulations. A single numerical run consists 
of generating a sequence of random flips such that the 
number of flips within a given time interval t follows a 
Poisson distribution, P„(t) = {'jt)^e~'^* /n\. The time 
dependent density matrix is then exactly solved by nu- 
merically integrating the von Neumann equation in small 
steps of At using the time evolution operators U and . 
The final numerical result is then obtained after produc- 
ing thousands of runs (5000 per dissipative curve in this 
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paper) each with a different RTN sequence, and then av- 
eraging density matrix over all sequences. This allows 
us to numerically simulate the quasi-Hamiltonian results 
which are inherently averaged over all RTN sequences. 
The results from the numerical simulation and quasi- 
Hamiltonian method are in exact agreement to within 
round-off error. 



IV. ANALYTICAL RESULTS FOR SHORT 
TIME EXPANSION AT ANY WORKING POINT 



Often tolerance levels for decoherence levels are quite 
stringent: a fidelity loss of more than one part in 10^ can 
destroy the result of a computation [62, £3^ . Hence e- 
folding times such as Ti and also the above Tenv are not 
necessarily the most relevant:it is important to look at 
short time behavior. We now analytically show that at 
short times, J always suppress decoherence for the Bell 
states {4'±y ^±) for an arbitrary working point 0, and for 
any set of noise parameters. 

The von Neumann equation can be expanded as follows 



(18) 



p{t) - p{0)-t / [H{t'),p{t')]dt' 
Jo 

= p(0)-i / [Hih),p{0)]dti 
Jo 

/ [H{t2)[H{h),p{0)]]dhdh 
10 Jo 



The ensemble average of p{t) has to be then calculated 
over all noise sequences 



A. Short time decoherence, for ^± at any 9 

The ensemble averaged non-zero Bloch vector com- 
ponents for the ip± Bell states, obtained using p{0) = 
\ip±){'<P±\ in the short time expansion of Eq. [l9j are 
listed in the Appendix along with further details. As 
seen in these equations, the Bloch vector components 
nox, n-xo, rioy and Uyo depend on J and do not contribute 
to the suppression of decoherence at pure relaxation 
{9 — tt/2) or pure dephasing {6 = 0) points. Whereas 
n^x and Uyy depend on J^. They are equal for pure 
dephasing and contribute to the suppression of decoher- 
ence at all working points except at pure relaxation. The 
Bloch vector components rizx, n^z, rizy and Uyz disappear 
at either pure relaxation or dephasing points. These cross 
terms along with n^y and Uyx do not contribute to the 
suppression of decoherence. 

The following approximate expression (for short times) 
can be obtained for the magnitude of the 'ijj± Bloch vec- 
tor by taking the sum of the squares of the above Bloch 
vector components 



3 - 8.g: 
2 



(21) 



where, = ^ ■ n'j was truncated after after the term. 
Note that the lowest order contribution from J is ob- 
tained only when Eq. 19 is expanded upto fourth order 



in time. From the functional form of Eq |53[ it is clearly 
seen that the decoherence will be suppressed with the 
onset of the exchange interaction. 



(Pit)) = 



N 



Jo 



{CN)dtidt2.-dtN (19) 



B. Short time decoherence, for <l)± at any 9 



where Cat is the N'^h order nested commutator, C^r — 
[H{tN), ■■■[H{t2), [H{ti),p{0)]]]. The ensemble average of 
(Cat) has to be then calculated over all noise sequences. 
This is done using the noise autocorrelation functions 
which are given by 



{s{t,)s{t2)-s{tN)) = {f\l[s{t,)^W,\t) (20) 



2 



_ g-27|t2fc-t2fc-l| 



fc=l 



where s{tj) = 1 for the first row of W and is —1 for the 
second row. Wj is as in EqjTjwith t — > tj+i — tj. Includ- 
ing a factor of {g/2)^ in the above expression accounts 
for the noise amplitude. 



Similarly, in case of the (j)± Bell state, p{0) = \4i±){4>±\ 
is used in the short time expansion of Eq. [19) As in 
the previous case, the lowest order contribution from J 
is only obtained when Eq |19| is expanded upto fourth or- 
der in time. The non-zero Bloch vector components are 
listed in the Appendix. As seen in these equations, the 
Bloch vector components nox,nxoTTLoz and Uzo depend 
on J, however at zero magnetic field none of these com- 
ponents will not contribute to the suppression of deco- 
herence in the short time expansion, rioz and Uzo respec- 
tively contribute only at pure dephasing and relaxation 
points. Whereas Uox and Uxo contribute only at interme- 
diate working points. Note that even though all of these 
components depend on the sign of J, the overall magni- 
tude of the Bloch vector (and hence the decoherence) is 
independent of the sign of J (just as in for n^{t)). 

In case of the 4>± Bell state, the following expression 
for the square magnitude of the Bloch vector in the short 



7 



time expansion of Eq. 19 



1 



2(Q 2 I 2 



(22) 
Bl)] 



where, = 
is seen in Eq 



n? was also truncated after term. 



It 



23 



that unhke the case of n^, (<) , J does not 
contribute in suppressing the decoherence at immediate 
times. Instead aids in maintaining the quantum co- 
herence in its place. Even though it appears that the 
J dependent terms will not suppress the initial decoher- 
ence due to their higher order time dependencies, our 
calculations show that the initial t^-decoherence can be 
compensated for with a sufficiently strong J at slightly 
longer times. 



QUASI HAMILTONIAN RESULTS FOR TWO 
FLUCTUATORS 



In this section we will demonstrate that even in the 
presence of two uncorrelated fluctuators coupled to the 
qubits, the exchange interaction still suppresses the de- 
coherence. 

The Hamiltonian for two qubits coupled to two uncor- 
related RTSs is, 

H=[^ + si(t)gi] ■ Si + [B + S2(t)g2] ■ Sa + JSi ■ Sa (23) 

where, gu and Sk{t) is the respective noise strength and 
RTN sequence for the fc*^ fluctuator. 6*^, is angle between 
gk and B and is the working point of the fc**^ qubit. The 
magnetic field is taken to be in the z-direction. 

Here we will only obtain analytical results for the case 
of pure dephasing for t/i^ . The two qubit dynamics is then 
once again reduced to an effective single qubit problem 
governed by = 2Jcr^ + [5ziSi(t) - gz2S2{t)](7x- The 
flipping probability matrix for two uncorrelated and un- 
biased fluctuators is W = Wi (g) W2, where Wi(2) has 
the same functional form as EqjTjbut with 7 replaced by 
7i(2). Analogous to Eq|9j the combined transfer matrix 
for the two fluctuators and the effective qubit at the j^^ 
instance of time is 



r. = w T 



(24) 



where T is a square matrix whose each column consists of 

the transfer matrices [T++,T_| ,T |_,T ]. Here, it is 

implied that T±,± = Tsj=±i,s2=±i = Lo + [JLz±gziLx± 
gz2Lx]At. Following the procedure outlined at the end 
of seclnlwe obtain the following quasi Hamiltonian 



(25) 



where, 

Hqa = i'yi{(^x - (^0) <8l (To (gl Lo + igzlUz <8l (To (gl La: (26) 



qb 



(To) <8 Lo + igz2(To ®(Tz®Lx (27) 

(28) 



iJ^Q, iJ^^ and i?^^ commute with each other. Hence in 
the absence of J, the time dependent Bloch vector can be 
solved for exactly. Analogous to Eq[8j for two unbiased 
and uncorrelated fluctuators, the time dependent Bloch 
vector is 



n(t) = ((/2|0(/i|)exp(-ziJ,i)(|n)® N2))n(0) (29) 

where, |ii(2)) — |/i(2)) — [ljl]/\/2 are the initial and 
final state vectors of the two unbiased fluctuators. This 
results in the following non-zero component of the Bloch 
vector 



2 r 



w - n 



k=l 



cos{ilkt) + — sin(17fct) 



p — ikt 



(30) 



where flk = \/9lk~ll- 

Next, for understanding the the effect of J on the Bloch 
vector, the matrix exponential in Eq|29|is calculated with 



the inclusion of H^^ using the Zassenhaus expansion (see 
Eq |12[ ). We derive the following approximate expression 
for the dephasing of the non-zero component of the effec- 
tive Bloch vector, valid only at short times or for small 
J or for small gzk 



<{t) « n 



Cfe cos{Jgzkt^) + sin(f2t) sm{2Jt) sm{Jgzkt^) 

(31) 



where Cfc = cos(51fct) -I- ^ sin(f2fct). As clearly seen from 
EqsjSO] and [STJ the Bloch vectors temporal dynamics in 
the presence of two uncorrelated fluctuators is simply a 
product of the single fluctuators states shown in EqsfTT] 
and [13] This suggests that the temporal dynamics of in- 
teracting qubits in the presence of multiple uncorrelated 
fluctuators (which can result in 1// noise) can be just as 
well understood in the single fluctuator picture. 

So far, in this section, we have analytically shown how 
J suppresses pure dephasing for ip± Bell state. The most 
general quasi Hamiltonian for two interacting qubits in 
the presence of two uncorrelated fluctuators with for any 
arbitrary set of working point, noise parameters, field 
strength and for any initial state is 



Hq — Hqy + Hqg + Hqb + HqJ 



(32) 



where 



Hq-/ = i['yi{(Tx — (To) <g (T0+ 12(T0 ® {(Tx — (To)] ® L'o ® L'o (33) 

Hqg = i[gi- L'] ® L'o®(Tz®(to~[- iL'o ® [gi ■ L'] (g) ao (8) 0-2(34) 
HqB = iBzL'o<S>[L',(SLo + L'o'S>L',] (35) 
Hgj = iJL'o(g)[L'^(g)Azx+Azx(E)L'^ + - (36) 

L'y ® Axz + Axz ®L'y+L'-,® Axy + Axy ® L'^] 
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Here, A^^; = (ctz 



+ cro) 

and Axy = ® (Tx - cry ay. L\^^^y^^ 
form of the S'0(3)-generators i^, whose first row and first 
column is padded with zeros. Lg is 4 x 4 identity matrix, 
9 = \9x,gy,g:^ and L' = [i^,L^,L^]. It is important 
to note that the quasi Hamiltonian 



is the 4x4 



Hq in Eq32 



64 X 64 matrix which when projected down using Eq 30 
acts on the 16 component Bloch vector instead of the 15 
component generalized Bloch vector. 

The decoherence and the disentanglement dynamics 
for any initial Bell state can be calculated for any given 
set of noise parameters, by exponentiating the above 
quasi-Hamiltonian (Eq 32 ) numerically. As an exam- 
ple, in FigjT] we have shown the dissipative dynamics of 
the V'+ entangled qubits in the Markovian, intermedi- 
ate and non-Markovian noise regimes for an equal mix- 
ture of dephasing and relaxational noise for both qubits 
i.e. 9i = 02 = 7r/4. Note that if the qubit working points 
are held equal, then the qualitative temporal behavior of 
dissipative process is the same as that of pure dephasing 
in all three noise regimes, but with overall longer decoher- 
ence times. As seen in the figure, the decoherence and the 
disentanglement is significantly delayed with increasing 
J in the non-Markovian and intermediate noise regimes. 
Similar to the single fluctuator case, in the Markovian 
limit the suppression of decoherence occurs when J 7. 
For the sake of simplicity, all noise parameters are held 
equal for both fluctuators in FigjT] However, J will sup- 
press the decoherence for any arbitrary choice of 9s, gs 
and 7s for both tj}± and (f>± initial Bell states. This can 
be easily verified numerically by using Eqs|32]and Eqj29] 



VI. INTERACTING QUBITS IN THE 
PRESENCE OF 1/f NOISE 

We now generalize our results for the interacting 
qubits, for a full spectrum of uncorrelated fluctuators. 
The general quasi Hamiltonian for two interacting qubits 
in the presence of n uncorrelated fluctuators, for any ar- 
bitrary set of parameters, with m {where m < n) fluctu- 
ators coupled to one qubit and n — ni fluctuators coupled 
to the second qubit is 



Hq — Hqy + Hqg + Hqb + Hqj 



where 



Ha 



Haa ~ 



L'o<»L'o 



HqB 
HqJ 



iY^{g~',-L')(8)L'o(g)ri''> + ... 

n 

iB^In ® (I/'z ® I/O + I/O ® L'j.) 

iJI„ ® {L'^ ® k^^ -f Az:^ ® I/L + ... 



(37) 

(38) 
(39) 



(40) 
(41) 



In is an identity matrix of dimension 2" and it is implied 
that 



'x(z} 



^0 



J") 



(42) 



The time dependent 16 component Bloch vector is now 
obtained from the projected temporal dynamics of the 
quantum and n classical TLSs as follows 



19^1=0.1, r^ = 0.1 ig^i=o.i,r^=i 




i 



50 100 150 
Time (1/B) 



50 100 150 
Time (1/B) 



50 100 150 
Time (1/B) 



FIG. 7: (color online) Bloch vector magnitude and 
concurrence as a function of time and Js for 1/)+ in the 
presence of two uncorrelated RTSs shown for (a) 
non-Markovian, (b) intermediate and (c) Markovian noise 

coupling regimes at 9i = n/4, and = 1. For all 
calculations shown in this figure, 51 = 52, 7i = 72 and 
01 =02. The the Bloch vector amplitudes are offset by 
+1 for clarity. 7, g and J are in units of B. 



n(i) = {fn\(g>-{f2\(S){fi\eM-iHqt)\ii)(g)\i2)...(S)\in)n{0). 

(43) 

As the dimension of the quasi Hamiltonian scales as 
2", calculating the coupled qubit dynamics in the pres- 
ence of a very large number of fluctuators can quickly 
become computationally very intensive. We therefore 
calculate the magnitude of the Bloch vector and the con- 
currence, using the quasi Hamiltonian in Eq |37[ for a 
set of eight fluctuators with 4 of them coupled to each 
qubit. A random 1// distribution of 7s is taken (rang- 
ing from 0.95 to 0.02), which is sufhcient for generating 
a 1/f noise power spectrum (where the power spectrum 
is S{uj) = J^^i/ili + '^^)) as shown in Figs. |8 (a) and 
(b). For the sake simplicity and to limit the parameter 
space we have held the noise strengths, gjS and their re- 
spective working point the same for all the fluctuators. 
The coupled qubits dissipative dynamics is shown in Fig. 
[Sjfor the ■0-1- state at 0i = 7r/4, for gj = 0.01, 0.1, 1. In 
all three cases it is seen that the exchange interaction, 
J, suppresses decoherence and disentanglement. As ex- 
pected the rate of decoherence rate itself increases with 
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FIG. 8: (color online) (a)Power spectrum, S{lj), for a 1// 
distribution of 8 uncorrelated fiuctuators. (b) Inset showing 
the noise power spectrum on a logarithmic scale. Magnitude 
of the Bloch vector and concurrence for the 7/)+ state shown 
as a function of time and J in the presence of the same 
fiuctuators with all giS held constant at (c) gj = 0.01 (d) 
gj = 0.1 (e) Qj = 1, at = 1. 9 = 7r/4 for all fiuctuators. 
The magnitude of the Bloch vector is offset by +1 for the 
sake of clarity and 7, g, J are in units of B. 



increasing noise strength, gj. For a full spectrum of fiuc- 
tuators, if gj is smaller than the smallest 7 then the noise 
due to the RTSs falls purely in the Markovian noise cou- 
pling regime, as shown in Fig. [8]-(c). If 5^ falls somewhere 
in between the selected range of 7s, then one has a mix- 
ture of Markovian and non-Markovian noise sources as 
shown in Fig. |8]-(d). However in this case, the Marko- 
vian noise sources tend to dominate and the oscillatory 
behavior, typically seen for non-Markovian noise, tends 
to get washed out. Whereas for the case of gj = 1 in Fig. 
[8]-(e), where one has a mixture of Markovian and inter- 
mediate noise sources, small oscillations superposed on 
top of a smoothly decaying function can be seen. This is 
similar to the dissipative behavior seen for a single qubit 
in the presence of broad spectrum noise |36] . 

As in the previous case of one and two fiuctuators, once 
the exchange interaction is turned on, the decoherence 
and disentanglement dynamics is suppressed for multi- 
ple fiuctuators as well. This suppression is proportional 
to the strength of J. However, for the 1// noise power 
spectrum, it is seen that the effect of J in suppressing the 
decoherence is apparent only when it is at least the same 
order of magnitude as gj . Finally, if all the uncorrelated 
fiuctuators are at the pure dephasing point, then the de- 
coherence and disentanglement dynamics of the coupled 



qubits are simply products of the dissipative dynamics 
due to each individual fluctuator. The short time behav- 
ior can be calculated for any number of fiuctuators with 
any arbitrary distribution of 7s and gs by taking the k 



sum in Eq 31 from 1 to n 



VII. SUMMARY 

In summary we have suggested a new way to prolong 
quantum entanglement and coherence via qubit-qubit in- 
teraction. We have analytically and numerically shown 
that the exchange interaction suppresses decoherence and 
disentanglement for the entangled Bell states. This is 
shown to be true for a single RTS, for two RTSs and for 
1// noise as long as the fiuctuators are uncorrelated. Our 
calculations are carried out using the quasi-Hamiltonian 
method, which is seen to be a particularly powerful 
method while dealing the many degrees of freedom asso- 
ciated with 1// noise. For the single fluctuator case, the 
suppression of decoherence is most apparent when J ^ 7 
and hence is more effective for non-Markovian noise. For 
1// noise, we conclude that J should be about the same 
order of magnitude as the strongest noise source, g, in 
order to alleviate decoherence. As the suppression of de- 
coherence at immediate times is key to performing high 
fidelity gate operations, we have analytically shown that 
the exchange interaction induced suppression of decoher- 
ence at short times is more effective for the ip± Bell state 
then it is for the 4>± state. This is true for any arbitrary 
qubit working point, fluctuator switching rate and noise 
strength. If however a large magnetic field is used, then 
the (f)± Bell state is more suitable for encoding quan- 
tum information. These results are vital for quantum 
information processing as they can be used to develop 
alternative methods (in addition to existing proposals) 
to enhance qubit lifetimes. 
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IX. APPENDIX 

In sec |IV[ for the 4'± Bell states, the ensemble aver- 
aged (pit)) is obtained by using p(0) = |V'±)(V''±I in the 
short time expansion (upto A^ = 4) of Eq. 19 For any 6, 



the lowest order contribution from J is obtained only 
when Eq. [19] is expanded upto fourth order in time. 
Note that the first order term {[H (ti) , p{0)]) = 0. Us- 
ing nij{t) = Tr [fji ® Uj X {p{t))], the following ensemble 
averaged non-zero Bloch vector components are obtained 
for 'ip± 
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_5 __10„ , ,4 

37t - 4 4 ^ ,„ ,3 

2 D T 2 ,4 

37t- 4 2 2^, „^,3 

y = -B^g^iyt - 2)t 



(44) 
(45) 
(46) 
(47) 
(48) 



^5^.5. [3t' - 27*' - (3Bf + 6j' + n')*-*] (49) 



g.g.l^t^ - It" + {B'i + 10 - J^')*-*] 



2 
3 

1 - U ~ ^74^ 



l-2(<?f + <?^)(t^-^7t^ 



-l + lgliit" ~2^t' -{Bl + ^f)t^] 



(50) 
(51) 

(52) 
(53) 



where, Vt = \J gi + gl - l"^ and = g^ tan(6'). 

Similarly for the (j)± Bell state, p(G) = |0±)(0±| is used 
in the short time expansion of Eq. [19] As in the case of 
V'i ; Eq[l9] has to be expanded at least upto = 4 to 



see the lowest order contributions from J at the pure 
relaxation/dephasing points. As expected, the ensem- 
ble averaged first order terms {[H (ti) , p{0)]) = 0. The 
resulting non-zero Bloch vector components are 



55! 



-2B^Jgxgzt'^ 
= 2B,Jgy 



(54) 
(55) 



[165," + gl{2- jt) + 4gl{3 - 2'yt)] (56) 

-^^B. [165," + 5^(8 - 5jt) + 45^(2 - 37*)] (57) 
8- 67* 



20 - 147* 
2 



8-67* 3 
5 B:,g^g^t 



-g^g, [3 - 27* - (3Bf + n^)t^] 

28 „2 

nxz 4- irg^g^B^ 



l-{2gl + %Bl)t'' ^--yglr" 



+ - [16B,* + (^gl + 24gi)Bl + gz^Q?] 



4 

2-37i 



l-2g^(*' + -7t + 



fB', + ln']t' 



2 .3 , B^,+ft\i 



t —n^ 
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